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Abstract. We compute the quark-gluon vertex in quenched lattice QCD, in the Landau gauge using an 
off-shell mean- field O (a) -improved fermion action. The Dirac- vector part of the vertex is computed for 
arbitrary kinematics. We find a substantial infrared enhancement of the interaction strength regardless of 
kinematics. 



1 Introduction 

In recent years, substantial progress has been made in our 
understanding of the nonperturbative correlation func- 
tions (propagators and vertices) of the fundamental fields 
of QCD and their relation to the phenomena of colour 
confinement and dynamical chiral symmetry breaking. It 
has emerged that at least in Landau gauge, a detailed 
knowledge of the structure of the quark-gluon vertex is 
essential for an understanding of the dynamics of quark 
confinement and chiral symmetry breaking as encoded in 
the quark Dyson-Schwinger equation (DSE), relating the 
quark propagator S(p) to the gluon propagator and the 
quark-gluon vertex /^(p, q), where p and q are quark and 
gluon momenta respectively: 



s-\ P ) = i^- 



m 



d 4 c 



(27T) 



n»S(p-q)D^(q)r u (p-q,q). (1) 



Here, D fIU (q) = P fIU (q)D(q 2 ) is the gluon propagator, with 
Piiv{q) the transverse projection operator. 

The overall shape of the gluon propagator is now quite 
well known, both from lattice QCD [HOIS] and from stud- 
ies of the coupled ghost-gluon Dyson-Schwinger equa- 
tions. It is now clear that if this is fed into the quark 
DSE together with a bare or QED-like vertex, the result- 
ing quark propagator will not exhibit a sufficient degree 
of chiral symmetry breaking [HEME] • 

The quark-gluon vertex is related to the ghost sector 
through the Slavnov-Taylor identity (STI), 



q^r^p, q) = G(q 2 ) (1 - B(q,p + q))S-\p) 



S-\p + q)(l-B(q,p + q)) 



(2) 



where G(q 2 ) is the ghost renormalisation function and 
B(q, k) is the ghost-quark scattering kernel. In particu- 
lar, if the ghost propagator is infrared enhanced, as both 



lattice [HE] and DSE studies [9] indicate, the vertex will 
also be so. This provides for a consistent picture of con- 
finement and chiral symmetry breaking at the level of the 
Green's functions of Landau-gauge QCD, where the same 
infrared enhancement that is responsible for confinement 
of gluons, provides the necessary interaction strength to 
give rise to dynamical chiral symmetry breaking in the 
quark sector. 

Confinement of quarks is still not fully understood in 
this picture, however. If we ignore the scattering kernel 
B in ([2|), dimensional analysis shows that the factor in 
square brackets is proportional to q as q — > 0, so 



r^q)~G(q 2 )~(q*y 



0. 



(3) 



Using the DSE result D(q 2 ) ~ (q 2 )~ 1+2K , the effective in- 
teraction in the infrared between a quark and an antiquark 
is 



V(q) = r(0,q)D(q)r(0,q) ~ (q 2 )~ 2K (q 2 )- 1+2 « = q~ 2 , 

(4) 

while a linearly confining interaction would be given by 
V ~ q~ 4 . This indicates that the quark-gluon vertex 
must contain an infrared enhancement over and above 
that contained in the ghost self-energy, connected with 
an enhancement of the scattering kernel. A recent DSE 
calculation [10] suggests that this may indeed be the case, 
with the running coupling a qg defined from this vertex 
diverging quadratically in the infrared. 

In two previous papers [TT| [T2] the quark-gluon vertex 
was computed on the lattice in two specific kinematics, the 
soft gluon ("asymmetric") point q = 0, and the quark re- 
flection ("symmetric") point q = —2p. Those results indi- 
cate a highly nontrivial infrared tensor structure, a result 
that has been qualitatively supported by semiperturbative 
DSE-based calculations [7|fT3j . However, these kinematics 
may not dominate the DSE (pQ). It is therefore necessary 
to compute the vertex in arbitrary kinematics, as far as 
possible. This is the focus of the present paper, although 
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at this point we will be restricting ourselves to the dom- 
inant, vector part of the vertex. Preliminary results were 
reported in [T4j . 

The structure of this article is as follows: In sec. [2] we 
briefly recap the formalism used in these studies. The main 
results are contained in sec. [3l while in sec. [4] we discuss 
the implications and possible future directions. Some tree- 
level lattice formulae are contained in the Appendix. 



2 Formalism 

We denote the outgoing quark momentum p and the out- 
going gluon momentum q. The incoming quark momen- 
tum is k = p + q. In the continuum, the quark-gluon 
vertex can be decomposed into four components Li con- 
tributing to the Slavnov-Taylor identity and eight purely 
transverse components Tf. 



r^p, q) 



4 

£ 

i=l 



\i(p 2 ,q 2 ,k 2 )Li^(p,q) 



(5) 



2=1 



The components Li and Ti are given in [TT] : 



L2,ll 


=- pp» , 


L3,n 


= -iP lt , 


La,^ 










=- n» , 




= fa? - qiv , 




— - i[q 2 cr^Pv 


^5, /Li 


— - icr^qv , 







(6) 



T 7 ,^ = - -(qP)a^ u P u - iP^a^xPukx , 



T 8 ,^ = - j^(Ji;\Pvk\- 



x+ ftp? 



(7) 



where P M = p^ + fc M , ^ = (pg)fc^ - (kq)p fI . In Landau 
gauge, for g 7^ 0, we are only able to compute the trans- 
verse projection of the vertex, I^f (p, q) = P^ v {q)r y {p^q)^ 
where P fJLU (q) = 5^ — q^q^/q 2 is the transverse projector. 
Since the vertex will always be coupled to a gluon prop- 
agator which contains the same projector, this is also the 
only combination that appears in any applications. The 
four functions are projected onto the transverse T^, 
giving rise to modified form factors 



Ai = Ai - q 2 r 3 , 

J2 



Ao — Ao 



-72 , 



(8) 



Ai 



Aa-yn, 



A4 = A 4 + q 2 r 4z . 



The lattice tensor structure is more complex, and ([5]) is 
only recovered in the continuum. The form factors also 



receive large contributions from lattice artefacts at tree 
level, so tree-level correction is required, as described in 
the Appendix. 



3 Results 

We have analysed 495 configurations on a 16 3 x 48 lattice 
at /3 = 6.0, using a mean-field improved SW action with 
a quark mass m ~ 115 MeV. This is part of the UKQCD 
data set described in [15] ; further details can also be found 
in [11]. 

The general lattice tensor structure, even for the Dirac- 
vector part of the vertex alone, is quite complicated, as 
the tree-level expression in eq. (fT6j) of the Appendix gives 
an indication of. This makes a determination of the full 
tensor structure of the vertex intractable with this lat- 
tice action. However, in the special case where both the 
quark and gluon momentum vectors are chosen to be 'per- 
pendicular' to the vertex, i.e. if we compute r fI (p^ q) with 

— % — 0? this structure simplifies considerably. There 
is no loss of generality provided rotational symmetry is 
restored in the continuum. We will here only study the 
leading, vector part of the vertex, as this is expected to 
have the cleanest signal, and tree-level lattice artefacts can 
be corrected multiplicatively. In continuum notation, we 
compute 



2 

Itr 7 ,r>, g) =(l-^)A' 1 



-[k 2 -p 2 -(kl-pl)]r 6 

■ K-i^-p^r^x". 



(9) 



(10) 



Of particular interest is the quark- symmetric limit, where 
the two quark momenta are equal in magnitude, p 2 = k 2 . 
In this case, tq is also eliminated, i.e. X{(p 2 ,q 2 ,p 2 ) = 
\[(p 2 ,q 2 ,p 2 )- Note that both the soft gluon (q = 0,p = k) 
and the quark reflection (q — —2p,p = — k) kinematics 
discussed previously in [TT|[T2] are specific instances of 
this more general case. In fig. [U we show X[ as a func- 
tion of the two remaining independent momentum invari- 
ants. The data become quite noisy as q in particular is 
increased, and also exhibit some 'spikes' and 'troughs' 
which at present we assume to be numerical noise and 
lattice artefacts. There appears to be some difference in 
the behaviour as a function of gluon momentum q and 
quark momentum p, but in view of the noise this cannot 
be further quantified. 

By interpolating thepoints in fig. [TJ we may reach the 
totally symmetric pointQ where p 2 = k 2 = q 2 . This kine- 
matics has a history of being used to define a momentum 
subtraction (MOM) scheme [16]. We show our results in 
fig. [2j together with the data from the soft gluon point 



1 As we are using antiperiodic boundary conditions in time 
for the fermions and periodic for the gauge fields, it is not 
possible to find a totally symmetric set of lattice momenta. 
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Fig. 1. The unrenormalised form factor A^ in the quark- 
symmetric kinematics p 2 = fc 2 , as a function of quark momen- 
tum p (long axis) and gluon momentum q (short axis). The 
form factor is given by the upper surface with values on the 
vertical axis ranging from 0.63 (blue) to 2.64 (red). The asso- 
ciated statistical uncertainties are given by the lower surface, 
with values ranging from 0.04 (blue) to 0.44 (red). 




P (GeV) 



Fig. 2. The unrenormalised form factor A^ at the totally sym- 
metric kinematics p 2 = k 2 = q 2 , as a function of the momen- 
tum p. Also shown is the same form factor in the soft gluon 
kinematics q 2 = 0. 



presented in [THIT2]. Again we find a strong infrared en- 
hancement, comparable to that of the two other kinemat- 
ics (see for example figs 7 and 8 in [12]). The qualita- 
tive behaviour of X[ as a function of p 2 is similar in all 
three cases, with a tendency to drop more rapidly as q 2 
increases, as one would expect if the vertex is infrared 
enhanced in all momentum variables. It does not appear 
possible to describe the three momentum directions with a 
single common variable, i.e. \[(p 2 ,q 2 ,p 2 ) = \i(t 2 (p 2 , q 2 )). 
Note that both sets should approach the same limit as 
p 2 0. 

Finally, figs. [3] and [4] show A" in general kinematics, 
for four different fixed values of q, as a function of the two 
quark momenta p and k. We expect all form factors to be 
symmetric in p 2 and k 2 (tq on its own is antisymmetric, 
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Fig. 3. The unrenormalised form factor A" for gluon momen- 
tum q = 0.555 GeV (top) and q = 0.877 GeV (bottom), as 
a function of quark momenta p and k. The values, shown by 
the upper surface, range from 1.0 (blue) to 2.2 (red), while the 
statistical uncertainties, illustrated by the lower surface, range 
from 0.03 (blue) to 0.15 (red). 

but is multiplied by p 2 — /c 2 ), and this is also what the 
figures show, within errors. The broadening of the data 
surface as q grows is simply a reflection of the increase in 
available phase space. 

The same qualitiative features as were found in the 
more specific kinematics, are reproduced here. At low q, 
we see a clear infrared enhancement, which disappears as 
q grows, reflecting the fact that at high momentum scales, 
only the logarithmic behaviour (which is too weak to be 
seen in these data) remains. At the same time, the level of 
the surface sinks, which reflects the infrared enhancement 
of A" also as a function of gluon momentum. 



4 Discussion and outlook 

We have determined the leading component of the quark- 
gluon vertex as a function of all three momentum invari- 
ants p 2 ,k 2 ,q 2 . We find an infrared enhancement in all 
momentum directions, although the quality of the cur- 
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Fig. 4. As fig. [3l but for gluon momentum q = 1.193 GeV 
and 2.321 GeV. For q = 1.193 the values range from 1.0 (blue) 
to 1.8 (red) while the uncertainties range from 0.03 (blue) to 
0.23 (red). For q = 1.193 the values range from 0.4 (blue) to 
1.5 (red) while the uncertainties range from 0.03 (blue) to 0.39 
(red). 



rent data does not make a further quantification of this 
feasible. 

These results have been obtained on a rather small lat- 
tice, and with a discretisation that gives rise to quite large 
tree- level lattice artefacts which must be corrected for. We 
therefore expect systematic errors to be quite large. To ob- 
tain more reliable results, and to extend this study to the 
full vertex structure at all kinematics, it would be desir- 
able to employ an action which is known to have smaller 
and more tractable tree- level artefacts. The Asqtad action 
has been employed successfully in computing the quark 
propagator [17], and unlike the SW action, only Ai and 
possibly A2 are non-zero at tree level, so no tree-level cor- 
rection will be needed for the remaining form factors. This 
action is also computationally relatively cheap, making 
large lattice volumes feasible. Another possibility is using 
overlap fermions, which have the advantage of retaining 
an exact chiral symmetry, which also protects all the odd 
Dirac components of the vertex at tree level. 
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Appendix 

We use the same notation as in [12], defining the lattice 
momentum variables 

K^p) = -sin(p M a); (11) 

Q M (p) = - a sin(p M a/2) ; (12) 

Cfj,(p) = cos(p M a) . (13) 
The tree-level vertex can be written 

r^(p, q) = c m [i F(p)A v (p) + B v (p)] x 

x { l^C^s) -iK^s) ^cr^K^j x 

i JK(k)A v (k) + B v (kj\ /D I (p)D I (k) , (14) 

where k = p + q;s = P/2 = (p + fc)/2; c m = 1 + b q ma; 
and Ay , By and Dj are defined in [12] . The transverse 
projection is given by 

as) 



Concentrating on the vector part of the vertex, we find 
that 

D I (p)Di( k) 



tr 



(7,/2 0) (P,9)) 



Q 2 M 

Q 2 (q) 



x A v {p)A v {k)K{p)-K{k) + B v {p)B v {k) 
2A v (p)A v (k)C fl (s)K fl (p)K^(k) 

x Y,C v (s)Q v {q) [K^p)K„(k) + K^k)K„(p) 

V 

A v (p)B v {k)K^s)K^p) + Avi^Bv^K^K^k) 
K{s) ■ Q{q) 



Hi) 



A v (p)B v (k)K fl (p) 
+ A v (k)B v (p)K^k)]Q^q) 



^C,q)lA v {p)B v {k)[K{p)-K{q)-K,{p)K ll {q) 



A v (k)B v (p) K(k)-K(q)-K ll (k)K ll (q) 



(16) 
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If we now choose to impose the condition — = 0, 
(fl6|) simplifies to 



IT 



D I (p)D I (k) 

B v (p)B v (k) 

(17) 



x | A v (p)A v (k)K(p)-K(k) + B v (p)B v (k) 



A v {p)B v {k)K{p)-K{q) 
-A v (k)B v (p)K(k)-K(q) 



The tree-level correction is carried out by dividing the 
quantity obtained nonperturbatively from the lattice, cor- 
responding to the lhs of (fTT|). by the expression on the rhs 

of CCD. 
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